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Abstract 

We present simulation results for long {N < 4000) self-avoiding walks in four dimen- 
sions. We find definite indications of logarithmic corrections, but the data are poorly 
described by the asymptotically leading terms. Detailed comparisons are presented with 
renormalization group flow equations derived in direct renormalization and with results 
of a field theoretic calculation. 
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1 Introduction 



Self avoiding walks (SAW's) are of practical importance since they form a model 
for randomly coiled linear polymers. But even more important for theoretical 
physics is that they represent in some sense the simplest critical phenomenon. 
More precisely, they formally are described by the m — limit of the 0{m) 
Landau-Ginzburg field theory Other members of this family are the Ising 
(m = 1) and the Heisenberg (m = 3) models. 

The most profound theoretical understanding of these critical phenomena is ob- 
tained by the field theoretic renormalization group evaluated near 4 dimensions. 
Above d = 4, all 0{m) models show mean field behavior. Below d = 4, the devia- 
tions of the critical exponents from their mean field values are of order e = 4 — d. 
This follows from the fact that the renormalized coupling constant is of order e 
in the infrared limit. The most extensively studied method for estimating critical 
exponents (the "e-expansion" ) involves resumming perturbation expansions in e. 
Precisely ai d = 4, the deviation from mean field behavior is given by logarith- 
mic corrections which can be predicted unambiguously by renormalization theory. 
Specifically, the leading behavior in the limit of long chains is ^ 

Rl ~ AT [log AT]", a = 1/4 (1) 

for the average squared end-to-end distance, and 

Cn ~ ^^[logAr]^ 13 = 1/4 (2) 

for the number of distinct walks. These predictions are basic results. In partic- 
ular, they do not suffer from any ambiguities inherent in a resummation of the 
e-expansion. It thus would be extremely useful if one could verify them by in- 
dependent means. The most obvious candidates for alternative calculations are 
exact enumerations of short walks and Monte Carlo simulations. Both methods 
have been applied in the past. 

Enumerations of chains with length up to A^ = 18 have been used in ^, ^ to 
verify eq.(2) {Rn was not computed in these papers), and recently chains were 
enumerated with A^ up to A^ = 21 ||^. From this the authors claim excellent 
agreement with eq.(2): the power of the logarithm in the best fit is /5 = 0.250 ± 
0.005. 

Monte Carlo Simulations, on the other hand, have been used previously for es- 
timating a in p. In addition, different exponents related to a and (3 were 
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measured in 0, If. While refs. 0, 0, |[ claim good agreement with eq.(l) (with 
e.g. a = 0.25 ±0.02 in [0) and with eq.(2), serious disagreement was found in 0: 
the best fit was obtained there with a = 0.31. Although no error bars were given 
in 0, the author obviously considered the value a = 1/4 to be ruled out. This 
seems a serious problem since the simulations of ref. ^ are by far more significant 
statistically than those of [§, 0, H, both concerning the chain lengths (A^ up to 
2400) and the sample size. Also, the method of analysis used in was justly 
criticized in since it introduces an uncontrolled bias and does not use the data 
optimally. 

The methods used in [0, |^ (in both papers, identical methods were used) are not 
easily compared to the present one. In different runs, these authors measured the 
average length {Nji)p of chains with fixed end-to-end distance R in grand canonical 
ensembles at p < pc {p is the fugacity), and the probability Q{j)) that two chains 
starting at neighboring sites never cross each other. We have not measured Q{p) 
ourselves (such measurements will be presented in [|10|); but we have measured 
{Nji)p. From this we shall argue in sec. 5 that the analyses of 0, |^ have large 
systematic errors, and that indeed deviations from mean field behavior is much 
larger than claimed there. 

In view of this situation (and since simulations of theta polymers in 3 dimensions 
[PI also gave discrepancies with the logarithmic corrections expected there 0]) 
we decided to perform simulations with much higher accuracy than those done 
previously. 

We find that the numerical results of |P are correct (though we do not agree 
with the conclusions drawn from them). In contrast, it seems that the findings of 
[^, 0, ^ were not completely uneffected by wishful thinking. This means that the 
leading logarithmic terms of (|^) and (^ are not sufficient to describe the behavior 
at any chain lengths which can be simulated in the foreseeable future (unless we 
assume, as was done in that the leading terms show a different power of log A^). 

The first corrections to the asymptotic laws (||) and (D) have been calculated by 
Duplantier |12[, but it turns out that even the corrected expressions are not fully 
consistent with the numerical data. These expressions involve a non-universal pa- 
rameter (an integration constant from integrating the renormalization group fiow) 
which should be the same for Rj^ and Cn- But for a good fit two different values 
have to be chosen for Rf^ and Cn, which is an internally inconsistent procedure. 
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However, from our data we can extract logarithmic derivatives dlnRf^/dlnN, 
dlnCN/dlnN, which may be compared directly to renormalization group flow 
equations of an appropriate direct renormalization scheme. Such a comparison 
does not involve unknown fit parameters. It thus is very pleasing that the flow 
equations consistently can be put into a form which compares very well to our 
Monte Carlo data. 

Stimulated by that success we use the original data to test another theoretical 
approach. As has been stressed previously |jT3|, |14|, the renormalization approach 
can be viewed as proceeding in two essentially independent steps. We first deter- 
mine a mapping from the physical model to its renormalized counterpart, and in 
a second step we use renormalized perturbation theory to determine the observ- 
ables, working completely within the renormalized model. It has been suggested 
that the mapping must be determined most precisely, but in the second step we 
may be content with low order perturbation theory. We find that this method 
allows for a consistent determination of the nonuniversal parameters and yields a 
good quantitative fit to our data. 

The paper is organized as follows: In the next section we present our Monte Carlo 
algorithm and the resulting data. Comparison with renormalization group flow 
equations derived in the spirit of direct renormalization is given in sec. 3, while 
comparison with field theoretic renormalization is presented in sec. 4. We conclude 
with a discussion of our results in sec. 5. 



2 Simulations 

Aiming at estimates for both Rj^ and Ctv, we decided not to use the pivot algo- 
rithm ||16|, though it should be the most efficient algorithm for estimating Rj^i in 



the limit — > oo. Instead, we essentially used a recursive and randomized im- 
plementation of the old enrichment method . This method is an improvement 



over the incomplete enumeration or 'recursive sampling' method used in |I8|, [19 
and was applied successfully in |jTl|, to polymer adsorption on surfaces, to theta 
polymers and to off-lattice polymers interacting via Lennard- Jones potentials. 

The basic structural element of incomplete enumeration is a subroutine STEP(x) 
which marks the site x as occupied and calls itself at all neighboring sites x ± ej, 
provided these sites are still free and provided A^ < N^ax- Before leaving the 
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subroutine, the site x is marked as free again. If each free neighboring site were 
visited with probabihty 1, this would give exact enumeration. In order to obtain 
a grand canonical distribution of walks where 

un = const CnP^ (3) 

is the average number of A^-step walks, one has to use a random number generator 
so that only a fraction p of all free neighbors are visited. 

The basis of the present improvement is the observation that it is sufficient to 
visit only one of the neighbors (provided it is free; otherwise the subroutine is left 
immediately), but in average (A/" — l)p times. Here A/" is the coordination number 
of the lattice, and we have assumed that we do not attempt any back steps as the 
corresponding site would not be free anyhow. The average chain length diverges 
when p Pc = Since (A/" — l)pc > 1 for all lattices, this means that we make 
> 1 attempts to continue each successful path if p ^ p^. Thus we choose a random 
neighbor (different from that we had come from), call STEP at this neighbor, and 
after having returned we call STEP again at the same neighbor with probability 



(A/" — l)pc — 1- The main difference to the algorithm of e.g., is that we always 
try at least once to continue the walk. 



Like the incomplete enumeration method of [18, 19 1 and like the method of [21 
the present method corresponds essentially to a random walk in the chain length 
with reflecting boundary conditions at = and N = N^ax- At Pc, it takes 
thus roughly N'^/D steps to obtain one statistically independent SAW of length 
A^. Here D is a diffusion constant which is of order unity in the algorithms of 



1^, |2T], |T9|. The main advantage of the present algorithm (apart from slightly 



shorter programs) is that here 

6={{Af-l)p,-l). (4) 

Since 6 ~ 0.033 for the simple hypercubic lattice in ci = 4, this gives a very large 
diffusion coefficient, and the method is more than one order of magnitude faster 



than those of [0, |2T|, |T9[. Eq.(^) can be understood as follows: in the present 
algorithm, a step back in is only needed when the walker hits an occupied site. 
The chance for this is ~ 6. In simulations near pc, this means that in average 
1/6 forward steps are made until one back jump reduces A^ by 1/5. In each of 
[O, El], [19], in contrast, the probability to make a back step is of order 1 (it is 



1 - pY for |T|, [g, and 1/(1 + pAf) for [|1|]). 
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We have simulated SAW's of length up to N^ax = 4000. Since it is not possible 
to store the sites occupied by such long walks in a simple bit map, we used 
a hashing procedure similar to (but somewhat simpler than) that used in []16| . 



We made only simulations very close to Pci where n^r is roughly independent of 
A^. As explained in [|1^, |T9[, this should be most efficient. Our total sample 
corresponded io tin 10^. This is also roughly the number of SAW's of maximal 
length Nmax, but all these walks are of course not independent. Often it is stated 
that this correlation between the walks is the main drawback of the enrichment 
method. In our version, it is not a big problem because of the ease with which 
walks are generated. The number of independent walks is given by the number 
of instances where the algorithm has reached A^ = N^ax between two returns to 
the main routine (corresponding to A^ = 0). Our sample contained ^ 1.2 x 10^ 
such independent walks. The total CPU time was ca. 800 h on a cluster of DEC 
ALPHA workstations. 

Our results for R%/N are shown in fig.l. In order to compare with eq.(l), we 
plotted its logarithm against In(lnA^). To demonstrate the importance of non- 
leading terms in logarithmic expressions, we show two slightly different quantities. 
For the solid line, A^ is defined as the number of bonds, while A^ is the number of 
sites for the dotted line. These two definitions differ by one unit and are clearly 
equivalent for N ^ oo. Nevertheless, we see that the difference is still important 
for A^ = 300! We show also the results of which have larger error bars but 
are otherwise in perfect agreement. From eq.(l) we expect our data to fall onto 
a straight line with slope 1/4. This is definitely not seen. Instead, the best 
straight line to the data (where A^ is interpreted as the number of bonds) has 
slope 0.311 ± 0.003. It would indeed give a perfectly acceptable fit. 

Our data for Cn are shown in fig. 2. More precisely, we there plotted CnP^ against 
In A^ for three different values of p. We also show the exact enumeration data of 
[^] with which we are in perfect agreement. But again we see hardly any sign 
of the predicted asymptotic behavior. We do not want to discuss in detail why 
several authors 0, ^, ^ were able to extract the correct asymptotic behavior from 
enumeration data (for noiseless data there exist very sophisticated methods to 
extract singularities), but obviously non- leading contributions are very large. 
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3 Higher Order Terms and Renormalization 
Group Flow 



Fortunately, the leading corrections to eqs.(|l],^ have been calculated in |12 

ajv = Rl/N = r[\n{N/a)f^ 



171n(41n(A^/a)) + 31 
^ 641n(A^/a) ^ "' 



Cj,/fi^ = c[ln(iV/a)]i/" 



171n(41n(A^/a)) -3 
^ 64 ln(A^/a) ^ ' 



(5) 



(6) 



The constant a has to be treated here as a free parameter. We thus have two 
parameters (r, a) to fit R]^/N, and two more (/i, c) if we also want to fit Cjy. 

The best fit to eq.(|) is obtained with r = 1.331, a = 0.1237. It is included in fig.l 
(dashed line). Over the range of interest it practically is undistinguishable from 
a straight line with slope 0.311. 

For Cn it is even more obvious that the leading term [ln(A^/a)]^/^ would give a 
very poor fit. A fit with eq.d), using a again as a free parameter, gives c = 
1.05, a = 0.055. This fit is shown as dotted curve in fig. 2. Notice that the values 
of a extracted from a^- and from Cjv differ considerably. No acceptable fit is found 
with a common value of a. 

Neglecting this problem for the moment, our data suggest that the critical value 
of p is Pc = 0.147622 ± 0.000001. This is to be compared to Pc = 0.147625 ± 
0.000002 as obtained in p. The values accepted in 0, || were substantially 
larger, 0.1490 ± 0.0003 and 0.1493 ± 0.0007, which partly explains why smaller 
logarithmic corrections were found in these papers (see sec. 5). 

We have thus been able to produce individual good fits to R%/N and to C^P^ , 
but the non-leading corrections are very important (in the latter case masking 
completely the leading behavior), and the fit parameters are not mutually con- 
sistent. This can also be seen by plotting the ratio {R%/N) / {C^p^) . Here the 
dominant terms cancel, and we obtain 

const [1 + ■ ■ •] (7) 



CnP^ 321n(A^/a) 

From fig.3 we see that this ratio is indeed a roughly linear function of 1/ ln(A^/a), 
provided we take a ^ 0.1. The slope of this function has the right sign but is 
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roughly twice as large as the value predicted by eq.(|^). To summarize, even in- 
cluding the leading correction terms we do not find a fully satisfactory explanation 
of the data. 

To proceed we note that eqs. are derived by integrating the renormalization 
group fiow equations, keeping only terms up to one loop order. It therefore is 
of interest to take a step back and compare the data against the more basic 
fiow equations. We first reconstruct the fiow equations from results given in the 
literature. 



Following [^], we start from dimensionally regularized perturbation theory in 
d = 4 — e dimensions. Denoting by b the bare coupling constant and by 



(27r)^/2 

a dimensionless coupling strength, we have up to second order in z |T5| 



N 



l + z{- + \)-z\l-^ + -) + 0{z' 
e I /e^ e 



(8) 



(9) 



(10) 



Obviously these expansions are singular at d — ^ 4. We thus introduce a renormal- 



ized coupling constant, which in the minimal subtraction scheme of reads 



8 2 .64 17, o 4, 



(11) 



From this we get the Wilson function |T^ 



Q 1 

W\zR, e] = N—zn\,>,, = -ezR - Azl + llz^ + . . . , 
in which we can take the limit e ^ without encountering any problems: 



W[zr] = -AzI + 17zI + 0{z_ 



R) 



(12) 



(13) 



While oi CAr/yU^ expressed in terms of zr still are singular, their derivatives 
with respect to InA^ are known to be finite. Indeed these derivatives yield fiow 
equations governing the renormalization of the chain length and of the partition 
function. Defining 



and 



0-1 N,e] - 



N 



N 



d\nN 



b,e 



b,t 



(14) 
(15) 
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we have [IT2l 

7 

ao[zR] = \imao[zR, e]= zr + -z\ + ... , (16) 
(Ti[zb\ = \miai[zn, e] = zr - 5zl + . . . . (17) 

Equations (^ and (|]) were obtained by first integrating eq.(|l2D at e = 0, yielding 
zji as a function of A^. This was inserted into eqs. ([T6|) resp. (p^Tj), and eqs.(|l^) 
resp. ([15|) were integrated again. During these manipulations, only the leading 
terms were kept, since higher order terms are not completely known anyhow. 

Avoiding these integrations we now directly compare our data to eqs. (p!2D -(|T7|). 
We first calculate the derivatives (p^ , (|15D as functions of A^. Of course, we have to 
replace derivatives with respect to by finite differences (we use A In = In 2 for 
first derivatives, and In 4 for second), but this should not present much problems 
in view of the slow variations of all functions involved. 



We thus show in fig. 4 
while 



.,iN.p)^±-l,,PjpE (19) 

In 2 Cn 

is plotted in fig. 5 for the same three values of p as before. From these plots, 
we should be able to obtain the same function zji{N) by inverting eqs.([l^),(p!^). 
A quick test shows that this is not so easy. The problem is that zr{N) turns 
out to be not very small (even for the largest values of A^), and the Taylor ex- 
pansions (p^,(p!7D obviously are poorly convergent. This is particularly true for 
(Ji, for which truncation in eq.(|17D after the quadratic term yields Ui < 0.05, in 
contradiction to the data for A^ < 30. 

A trick which helps - though its justification is far from obvious - is to change 
the expansions in eqs.(|T4|),(p!7|) into Fade approximants, 

Mzr] = -. — > = , ■ (20) 

1 - 7zr/2 1 + bzR 



These can be inverted to give 



- - Zr l^-LJ 



with 

1 + 7(To/2 ' ^ 1 - 5ai 



4"' = , 4^' = ^ . (22) 
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In fig. 6 we have plotted zj^"* and z^^^ as obtained from tlie finite-difference ap- 
proximations. We see still some disagreement for small which may be due to 
higher order terms in ctq [-Zij] , cxi [^ij] and/or to corrections. But for large N 
we find very acceptable agreement, provided we take p = l/^ = 1.476223. 

The final test of the theory consists in checking whether these estimates of zr 
satisfy the differential equation N dzji/dN = W[zr] or, rather, its finite-difference 
approximation. Again we find problems due to the slow convergence of the Taylor 
expansion for and again we take recourse to a Pade approximant, 

W[zr] = . (23) 

^ 1 + 17zr/4: ^ ' 

In fig. 7a we have plotted this form of W with argument z^^'' against N , together 
with 

A4"^(Ar) ^ 4")(2Ar)-4")(Ar/2) 
AlnAT ln4 ■ ^ ^ 

We see very nice agreement, in particular for large values where it is most 

significant. For small A^ agreement could indeed be improved by adding a constant 
of 0{1) to A^ so that a^v is replaced by Rj^/^N + 0.16), see fig. 7b. Such a 1/A^ 
correction would also improve the agreement in fig. 6 at intermediate values of A^. 

We thus conclude that our data are fully consistent with the renormalization 
group flow predicted theoretically, though there is some ambiguity related to the 
use of Pade approximants and - to a much lesser degree - to the treatment of 
1/A^ corrections. We want to stress that this analysis is particularly interesting 
since, except for the size of the nonuniversal 1/A^ corrections, it does involve no fit 
parameter! Note further, that this analysis uses only the information underlying 
also eqs.(|) or (|^). This suggests that solving the renormalization group equations 
to one loop order and substituting the result into expressions for higher order 
corrections we lose meaningful information which still can be identified on the 
level of the flow equations. 



4 Field Theoretic Renormalization 



As mentioned in the introduction, the renormalization program may be separated 
into two steps: We first map the physical 'bare' model on a renormalized model 
which exhibits the scale invariance of long chains. We then calculate Rj^, Cn 
perturbatively within the renormalized model. 
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The bare model depends on the couphng constant b, the chain length N, and a 
microscopic length scale / which is of the order of the lattice spacing in the com- 
puter experiments. In the renormalized theory these parameters are replaced by 
the renormalized coupling u, the renormalized chain length Nr, and a renormal- 
ized length scale Ir. The mapping takes the general form 

6 = const [ — ] Zu{u,1/Ir)u (25) 



2 

I 



iV= f^) Z^{uJ/Ir)Nr, (26) 



where the renormalization factors are normalized according to 

Z„(0, 1) = ZniO, 1) = 1 . (27) 

The observables of interest here can be expressed as 

R% = 2dllNRAR{u, Nr) (28) 
Cn . Z{u,1/Ir) 

~N = const . Ac{u, Nr) , 29 

/i^^ Zn{uJ/lR) 

where Z{. . .) denotes another renormalization factor and the amplitudes Ar, Ac 
are to be calculated by renormalized perturbation theory. In general the renor- 
malization factors have to be chosen to absorb the leading dependence on the 
microscopic scale I: for fixed Ir, all observables have to become independent of 
/ up to corrections of order t^/l\NR ~ Using the scheme of dimensional 

regularization, defined by taking the limit / for d < A, factors such as to 
make the renormalized theory finite in four dimensions. 

The length scale Ir is a free parameter of the renormalized theory. Under a change 
of Ir the other parameters change according to the flow equations 

eu - f3{u) (30) 
= Cn{u) (31) 

OlnlR-^^^^^ ^''^ 
where all derivatives have to be taken at fixed b, I and e. Based on Pade-Borel 
summation of higher order perturbation theory, a parametrization of the flow 
functions has been given in p4| : 

5 

Cn{u) = —u + -u^ — aiu^ + a2U^ (33) 



OIyiIr 

5 In Zn 
dlnlR 
dlnZ 
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C(m) = + Asm' (34) 

/5 M = Au^— (35) 

1 + a^u 

with the constants ai = 3.6328, 02 = 3.8953, 03 = 0.04395, 04 = 1.555 and 
as = 3.5962. 



(Our notation differs somewhat from that of [El]: u = 8^1^31, = ~C 



m , - M 



C = Cf^~^- We further note that recently some mistake has been found in 
the five loop contributions to the renormahzation factors. However, due to the 
remarkable stabihty of the Pade-Borel results this is not expected to yield serious 
changes in the parametrization of (^)-(l35|).) 

Starting from arbitrary initial conditions set aX Ir = I we now integrate the flow 
equations to find functions u{Ir) e.t.c. in analytic form. We then fix the final 
value of Ir by the condition 

Nr = 1, (36) 
which implies that Ir is of order R^. There result the expressions 

gl.951u-1.422u2+o.734u^ 

= ^%V4(i + i.555u)i-369^«^^' 1) ' ^37) 

g2.101u-1.434u2+0.734w3 

The renormalized coupling constant is implicitly determined by the equation 

AT = n, ^^ + 1^555n)'-''% jL,(i.951n-1.422n^+0.734n3) _ (gg^ 

Here /q, cq, no are nonuniversal fit parameters absorbing the initial conditions in 
the integration of the flow equations. 

So far we have been concerned with the renormahzation group mapping. Calcula- 
tion of the amplitudes is a straightforward exercise in renormalized perturbation 
theory. It yields 

11 

Ar{u,1) = 1-- + 0{u') (40) 
Ac{u,l) = l + - + 0{u') (41) 

Eqs.(p7|)-(|^) constitute our final result. From the derivation it should be clear 
that the calculation of the amplitudes is a problem well separated from the de- 
termination of the renormahzation group mapping. Indeed, within the minimal 
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subtraction scheme it is the singular terms in expressions like eq.(||),(|l^) which 
determine the mapping, whereas the amphtudes result from regular terms. 



We now compare to the Monte Carlo data. Fitting aj\j in the range 300 < N < 
4000 we find Iq = 1.0436 and = 0.275 (see fig. 8). The same value for uq together 
with Co = 0.864 would also give an acceptable fit to Cjy. But the best fit to the 
latter (see fig.9a) is obtained with no = 0.316, cq = 0.870, and 

= 0.1476223 ± 0.0000005 . (42) 

The latter value for no gives also an excellent fit to ajsj, provided is replaced 
by + ni in the definition of on, with rii = 0.56 (see fig.9b). In the latter 
fit, Iq = 1.0466. A priori, such a change would be well within the uncertainty 
intrinsic in the definition of A^. But fig. 7 and analogous plots of the dependence 
of u on Ifj in the present scheme suggest that such a value of rii is somewhat large. 
We must note, however, that introducing rii we can account for 1/A^ corrections 
only in a very crude way PB|. We thus find agreement with the field theoretic 



prediction, but there remains some room for further improvement in the region of 
short chains. 

It is of some interest to note the range of the renormalized coupling resulting from 
our fit: 0.067 < u < 0.15 for A^ in the range 4000 > A^ > 100. Thus even for 
A^ ^ 4000 u is not particularly small. On the other hand, it is sufficiently small 
for the parametrizations of p4| to be justified. 



5 Discussion 



We have seen that self-avoiding walks on the hypercubic lattice in four dimensions 
show clear logarithmic corrections. These corrections are in perfect agreement 
with the predictions of the renormalization group, though they would be very 
poorly described by the asymptotically leading approximation. Thus claims |^, 
^ that exactly these leading terms have been seen in simulations must be viewed 
with some reservation. On the other hand, our numerical results agree perfectly 
with simulations reported in |^, although the conclusion drawn in that paper 
seems to be wrong. Though our data are perfectly fitted by power laws in InA^ 
— in particular, 

R%/N ^[\nN]^-^\ (43) 
these do not represent the asymptotic behavior. 
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In order to understand the discrepancy with [0, ^, we have repeated their com- 
putation of {Nji)p and their subsequent analysis, but with several modifications: 

(i) using our algorithm, we calculated {Nji)p not only for a single fixed value of R 
{R = 7m 0, i), but for all R < ^500; 

(ii) due to the efficiency of our algorithm (and improved hardware) we have much 
higher statistics; 

(iii) we use data much closer to the critical point: Pc — P > 0.0012 compared to 
Pc — P > 0.018; and 

(iv) in the analysis we use our very accurate estimate for pc- 

Assuming that the decay of the two-point function is dominated by the smallest 
mass m in the model (the inverse correlation length ^), the authors of 0, |^ obtain 

{Nr)p ^R— (44) 

with 

m= l/e~ V^~^[ln(p,-p)]-^, AT =1/8. (45) 

Thus, {Nr)p/R should be a function of p alone, independent of R. Our data 
shown in fig. 10 clearly demonstrate that this is true only for R>> C,, while 

{NR)pr^RlnR (46) 

for i? < ^. This is not surprising: the decay of the correlation function should 
be dominated by the smallest mass only for R» C,, while effects of the effective 
coupling should be visible at shorter distances. Indeed, keeping R fixed while 
p Pc we encounter the region where the data should be analysed via a short 
distance expansion, and eq.(|^) is no longer justified. According to fig. 10, R = 7 
is not yet in the regime where eq . (^4]) holds for the range of p considered in [0. 
This gives already a first indication that this coupling is larger than anticipated 
in [0 (we might add that the fiuctuations seen in fig. 10 are not statistical, but are 
lattice effects). 



In fig. 11 we show our data for R = 7, vTSO and ySOO together with those of [0, 
and with the mean-field prediction 

On the one hand we see that the slopes deviate strongly from mean field prediction, 
much more than they differ between different values of R. On the other hand. 
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we see the dramatic effect of a wrong choice of pc- Indeed, while we use our 
value of Pc when plotting our own data, we show the data of twice: once 
plotted using our p^, and once using their own estimate of Pc- We see that the 
latter reduces the deviation from mean field considerably. Thus we see that the 
logarithmic corrections are indeed larger than estimated in (in spite of the fact 
that the too small value of R alone would have led to their overestimation!), but 
the systematic uncertainties and the lack of higher order predictions from field 
theory prevent a more detailed analysis. The same comments should hold for 
the data of as this author used the same values of R and of p as 0, and his 
estimated Pc was even further from our estimate. 

To what chain lengths would we have to go in order to see the asymptotic behavior 
Rl/N, CN/fi^ ~ [In AT] 1/4? Our renormalization group calculation can easily be 
extended to larger and show that eq. (|i3| ) no longer would give a good fit above 
= 10^, but even at = 10"^ the effective exponent is ~ 0.285 instead of 1/4. 

The analysis in the field theoretic framework stresses the importance of a good 
quantitative form of the renormalization group mapping. For the amplitudes then 
a low order approximation seems sufficient. This supports previous findings, both 



in polymer physics |13| and in the physics of liquid helium |T3]. If a good form 
of the mapping were missing, this would justify the often expressed pessimism 
about the possibility to see logarithmic corrections. Such a situation prevails 
e.g. for ^-polymers in three dimensions. There only leading terms in the sense 
of eq.(|^),(|^) have been computed and they disagree badly with simulations 
[0. Although the disagreement there is even worse than the disagreement with 
the leading terms in the present case, one might suspect that also there the poor 
representation of the renormalization group mapping employed is at the root of 
the problem. 

The present analysis showed large subleading corrections because of the rather 
large value of the renormalized coupling constant, even for our longest chains. 
This raises the question whether the asymptotic behavior could be seen for much 
shorter chains in models with weaker repulsion between chains. Two such models 
come immediately into mind: chains with attraction between neighboring sites 
{6 polymers above Tg), and the Domb- Joyce model [^. In the latter, n > 1 



monomers can occupy the same site, but the contribution to the partition function 
gets a weight (1 — w)""^ for each multiple occupancy. For sufficiently small w one 
is arbitrarily close to the free case, and also the renormalized coupling constant 
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should be small. 
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Figure Captions: 



Fig.l: Plot of the logarithm of the swelling ratio, ln(i?^/A^), against In(lnA^). 
The solid line are our data with taken as the number of bonds; the dotted 
line represents these data with as the number of sites. The diamonds are the 
data from and the dashed line is the fit with eq.@. The dashed-dotted line 
indicates the slope predicted by the leading term given in eq.(l). The statistical 
errors of our data are roughly oc They are thus largest for = 4000, where 
AlnR% = 0.0005. 

Fig. 2: Semi-logarithmic plot of C^p^ against InA^. The plot shows our MC 
data for 3 different values of p very near p^, the fit with eq.(^, and the exact 
enumeration data from [^]. The statistical error of our data for = 4000 is 
ACn/Cn = 0.0013. 

Fig.3: F\otof{R%/N)/{CNP^) against l/ln(10A^). According to eq. (0) a straight 
line with negative slope is expected for p = pc, with a slope as indicated by the 
dotted line. 

Fig. 4: Function o"o(A^) against A^ (full line). Also shown is z^^\n) as obtained 
by the Pade approximant eq.(|^) (dashed line). 

Fig. 5: Function cri(A^) against A^, for the same 3 values of p as in fig. 2. 

Fig. 6: Functions z]^\n) (full line) and z]^^ (dashed). The latter is again shown 
for the same 3 values of p. 

Fig. 7: Panel (a): function W^f^;^] (full line) and the "derivative" of with 
respect to In A^ (dashed line). Both curves should agree up to higher orders in z^. 
Panel (b) shows the same data, but in the definition of we have replaced A^ 
by A^ + 0.17. 

Fig. 8: Fit to the swelling ratio with Iq = 1.0436 and hq = 0.275. 

Fig.9: (a) Fit to with cq = 0.870 and uq = 0.316; (b) fit to with the 

same no = 0.316, with Iq = 1.0466, and with rii = 0.56. 

Fig. 10: Average chain lengths at fixed p against their end-to-end distance R. 
Actually, the ratio {Nji)p/R is plotted. The values of ln(l/p) (from top to bottom) 
are: 1.9143, 1.915, 1.916, 1.918, 1.921, 1.925, 1.932, 1.948, 1.966, 1.984, 2.002, 
2.021, 2.066,2.100, and 2.150. 
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Fig.ll: Full lines: {Nr)p/R against pc - p for R = 7, 12.25, 22.36 (from bottom 
to top); dotted line: slope {Nji)p oc (pc — p)^^^"^ predicted by mean field theory; 
diamonds: data of [0] (for R = 7), using our value of Pc] crosses data of using 
their value of pc- 
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